ANALYSIS OF A ONE-DIMENSIONAL PRESCRIBED MEAN CURVATURE 
EQUATION WITH SINGULAR NONLINEARITY 

NICHOLAS D. BRUBAKER* + AND JOHN A. PELESKO** 

Abstract. In this paper, the classical solution set (A, u) of the one-dimensional prescribed mean curvature 
equation 

- ( , " ... ) = T . X u > ~ L < X < L > u(-L) = u(L) = 0, (*) 

\y/l + (u') 2 J (1-U) 2 

for A > and L > 0, is analyzed via a time map. It is shown that the solution set depends on both parameters 
A and L and undergoes two bifurcations. The first is a standard saddle node bifurcation, which happens for all 
L at A = A*(L). The second is a splitting bifurcation, namely, there exists a value L* such that as L transitions 
from greater than or equal to L* to less than L* the upper branch of the bifurcation diagram of problem (*) 
splits into two parts. In contrast, the solution set of the semilinear version of problem (*) is independent of L 
and exhibits only a saddle node bifurcation. Therefore, as this analysis suggests, the splitting bifurcation is a 
byproduct of the mean curvature operator coupled with the singular forcing. 
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1. Introduction 

The study of nonparametric surfaces of prescribed mean curvature, i.e., the study of solutions 
of 

-div^= = Xf(u), ieO; u = 0, x e 90, M l) 

goes back to 1805 and 1806 when Thomas Young [64] and Pierre-Simon Laplace [42] separately 
looked at the properties of capillary surfaces, namely, where / is linear. Later the theory, which 
still focused on capillary surfaces, was put on a solid mathematical foundation by Gauss [30] and 
attracted the attention of many nineteenth century scientific luminaries. During the first half of 
the twentieth century, the problem fell out of vogue, taking a backseat to the study of constant 
mean curvature surfaces (c.f. [56]); however in the past few decades, due to mathematical 
advances and the miniaturization of technology, the study of problem (1.1) has come back into 
focus. Most modern authors have focused on the existence, nonexistence and multiplicity of 
positive solutions. The case where f2 C M n , for n > 2, has been studied by numerous authors: 
Concus and Finn [19, 20, 21, 22, 23, 24, 25]; Giusti [32, 33, 34, 35]; Gilbarg and Trudingcr [31]; 
Ni and Scrrin [49, 50, 51]; Finn [28, 29] (and the references therein); Peletier and Serrin [61]; 
Atkinson, Peletier and Serrin [4, 3, 5]; Scrrin [62]; Ishimura [38, 39]; Kusano and Swanson [41]; 
Nakao [48]; Noussair, Swanson and Jianfu [52]; Bidaut-Veron [10]; Clement, Manasevich and 
Mitidieri [17]; Coffman and Zicmcr [18]; Conti and Gazzola [26]; Amster and Mariani [1]; Habcts 
and Omari [36]; Le [43, 44]; Chang and Zhang [16]; del Pino and Guerra [27]; Moulton and 
Pelesko [46, 47]; Bcreanu, Jcbelean and Mawhin [8, 9]; Obersnel and Omari [54, 55]; Brubaker 
and Pelesko [14]; Brubaker and Lindsay [13]. Also, the case where n = 1 has been studied by 
numerous authors in a recent series of papers: Kusahara and Usami [40]; Bencvieri, do O and 
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de Medeiros [6, 7]; Bonheure, Habets, Obcrsnel and Omari [11, 12]; Habets and Omari [37]; 
Obersnel [53]; Pan [57]; Li and Liu [45]; Burns and Grinfeld [15]; Pan and Xing [59, 58]. One 
of the fascinating aspects of problem (1.1) that these studies have revealed is the disappearing 
solution behavior shown to be present for numerous choices of / (cf . [59] , [58] and the references 
therein). For example, Pan showed in [57] that the solution set of 

/ „.i \ 

= Ae", < x < L; u(0) = u(L) = 0, (1.2) 



which is a quasilincar analogue of the semilinear Gelfand-Bratu equation, is fully characterized 
by the following theorem for A > (see Figure 1 for the resulting bifurcation diagram). 

Theorem (Pan [57]). Assume that A > 0. 

(i) If L < 7T, there exists constants A* and A* such that (a) for any A G (0, A*) U {A*} there 
is a unique positive solution of (1.2); (b) for any A G [A*, A*), there exists exactly two 
positive solutions of (1.2); (c) for any A € (A*,oo), no positive solutions of (1.2) exist. 

(ii) If L > n, there exists a constant A* such that (a) for any A G (0, A*) there exists exactly 
two positive solutions of (1.2); (b) for A = A*, there is a unique positive solution of (1.2); 
(c) for any A G (A*,oo), no positive solutions of (1.2) exist. 




Figure 1. Bifurcation Diagrams of (1.2) for L < n and L > n. Note that in (b), the bifurcation diagram 
continues; in particular, A — > 0+ -> oo [57]. 

In this paper we continue to study the disappearing solution behavior of prescribed mean 
curvature problem (1.1), with n = 1, by considering the case where / is an inverse square 
nonlinearity that is singular at u — 1. Namely, we study the solution set of 

with u < 1 in [—L,L], for positive A and L. Previously, these types of equations (i.e., one- 
dimcnsional singular prescribed mean curvature equations) have been studied in a general 
context by Bonheure et al. in [11]. In particular, applying their results to problem (1.3) gives 
the following: 

• if A > is sufficiently small, then there exists at least two classical solutions of problem 
(1.3) [11, Theorems 3.1 and 5.1]; 



ANALYSIS OF A ONE-DIMENSIONAL PRESCRIBED MEAN CURVATURE EQUATION 



3 



• there exists a A* > such that for all A > A*, no (classical or non-classical) positive 
solutions of problem (1.3) exist [11, Theorem 6.1]. 
To build on their work, we use the method of time maps to give a description of the exact 
number of classical solutions, i.e., C 2 (— L, L) n C[— L, L] solutions, when the parameters A and 
L vary. In doing so, we prove the following theorem that fully characterizes the solution set of 
(1.3). 

Theorem 1.1. Let L > and L* be defined by L* := m&x\ >0 g(\), where 

ir3 ^(lo g 1± 4 ZI "-V / ^), /-A, (0,1), 
{ 1 :!. ' for\ = l, ( L4 ) 

( A 2_ A 1)3 /2 (v / A^T-sec- 1 (A)), /or A € (l,oo). 

(i) If L < L* , then there exists three values A*, A** and A*, which depend on L, such that 

(a) for A g (0, A*] U [A**, A*), (1.3) has exactly two positive solutions; 

(b) for A g (A*, A**) U {A*} ; (1.3) has exactly one positive solution; 

(c) for A > A*, (1.3) has no solutions. 

(ii) If L > L* , then there exists a value A*, which depend on L, such that 

(a) for A g (0, A*), (1.3) has exactly two positive solutions; 

(b) for A = A*, (1.3) has exactly one positive solution; 

(c) for A > A*, (1.3) has no solutions. 

Furthermore, in both cases, A* < min {L 1 , 7r2 /(27I/ 2 )} . 

Remark 1.2. Using (1.4), the value of L* can easily be approximated: L* w 0.3499676. 

The results of this theorem are illustrated in Figure 2 for the two separate cases. Conse- 
quently, we see that the solutions set of (1.3) contains two bifurcations: 

(1) a saddle node bifurcation occurs, for all L, at the point (X*(L), \\u(- ; A*, L)||oo), where 
u(x; A, L) is the solution of (1.3); 

(2) a splitting bifurcation occurs at L = L* . That is, when L > L* , the upper solution 
branch of the bifurcation diagram of (1.3) is continuous (see Figure 2(d)); however, 
when L < L* , the upper solution branch splits into two parts (see Figure 2(c)). 

Remark 1.3. Equation (1.3) arises in the study of microelectromechanical systems, namely, [14] 
proposed that the steady state deflection, u, of an elastic membrane due to an electrostatic 
forcing satisfies 

— div — ^ U = = — — — ttt , ieO: u < 1, x g fl; u = 0, x g dfl, 

^TT^WW (i - u) 

where < e <C 1 is an aspect ratio of the device, (i is a positive ratio of the reference electrostatic 
force over the reference elastic force and O is a open, connected, bounded domain in R™, for 
n = 1,2. Therefore, when n = 1 (without loss of generality, we assume that the domain is 
centered at the origin), the change of variable y — e~ x x yields (1.3), where A := e 2 /i and 
L:=l/e. 

Remark 1.4. The disappearing solutions behavior exhibited in case (i) of Theorem 1.1 is due 
to the mean curvature operator; specifically, from [60], we know that for the related semilinear 
problem 

- u" = -L<x<L; u{-L) = u(L)=0, (1.5) 

(l u) 
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(a) L < L* (b) L > L* 




(c) L < L* (d) L > L* 

Figure 2. (a) Surface where the points on it correspond to the solutions of (1.3). Here, the black contours 
represent solutions for fixed L, which yield bifurcation curves that capture the qualitative shape described in 
Theorem 1.1; namely, for L = 3/10 and L = 6/10, which coincide with case (i) and case (ii), respectively (see 
Subfigures (c) and (d)). (b) Contour plot of the surface given in Subfigure (b). The new lines — (thick, dashed), 
(thick) and (dashed) — are the locations o/A*(L), A**(L) and\*(L) of Theorem 1.1. (c) Bifurcation diagram 
of (1.3) for L = 3/10; (d) Bifurcation diagram of (1.3) for L = 6/10. 



there exists a A* > such that for A <G (0, A*) there are exactly two solutions of (1.5), for 
A = A* there is a unique solution of (1.5) and for A > A* there are no solutions of (1.5). So, in 
other words, the bifurcation diagram of (1.5) depends on only a single parameter, which can 
easily be seen through the change of variable y — ^fkx, and has same qualitative shape as the 
one shown in Figure 2(d). Hence the disappearance of solutions is a byproduct of the mean 
curvature operator coupled with the singularity. 

We start the next section by proving some necessary conditions for the existence of classical 
solutions of (1-3). From this, we restrict the parameter space to where solutions may exist and 
use the inherent symmetry of (1.3) to reduce it to an equivalent problem. Then using these 
restrictions, in Section 3, we introduce and analyze a time map a la [57] and [63]. From this 
analysis, we then prove Theorem 1.1. 
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2. Necessary conditions for the existence of solutions 

In this section, we derive some properties about the solutions of (1.3). Also, we show that 
solutions cannot exist under certain conditions. 

Lemma 2.1. If u(- ;X,L) is solution of (1.3) for fixed A > and L > 0, then u > in 
(-L, L); furthermore, u is strictly concave down in this region, which implies that its maximum, 
IMloo G (0)1)> * s unique. 

Proof. From carrying out the differentiation in differential equation of (1.3), wc obtain u"(x) < 
for all x in (—L,L) and A > 0, which gives the desired result. □ 

Next, it is easily verified that the ordinary differential equation given in (1.3) has the first 
integral 

■ 1 -t^ = E, (2.1) 
y/1 + (W) 2 1 ~ u 

where E is a conserved quantity of the system and can be determined by the maximum deflection 
of u. That is, let c in (-L, L) be the unique value such that u(c) = |M|oo; hence, it'(c) = and 
first integral (2.1) gives E = 1 — A(l — ||it||oo) _1 ) which implies 

, 1 -^- = 1 * (2.2) 

y/1 + {U'f l-U 1-IMI 



| OO 



Using this we can prove the following necessary condition for classical solutions. 

Theorem 2.2. Let u(-;X,L) be a solution of (1.3) for A > and L > 0. Then ||u||oo < 
1/(1 + A). 

Proof. For contradiction assume that ||w||oo > (1 + Hence, 

1 - A - ||u||oo 

From Lemma 2.1, we know that u attains its unique maximum at a value, say c, in (—L,L), 
i.e., u(c) = IMloo > u(x) for all x in (— L, c) U (c, L); therefore, inequality (2.3) implies that 
there exists an x in (— L, c) U (c, L) such that 

= i-a + ^)Ni~. 

1 - A - Hulloo 

Now, from first integral (2.2), we know 

'( \T\ 2 — (1 — m) 2 (1 — Halloo) 2 , „ 

u[x; ' j + [(i-«)(i -Hoc)- xiWul^-uW X ^ X - 

Since v! is continuous and bounded in (—L,L), we take x — > x so 

(l-U^l-IHloo) 2 



u'(x; \,L) 2 = -l + lim 



[(l- U )(l-|| U || 00 )-A(|| U || 00 - U )]2' 
however, the limit on the right-hand side diverges to infinity and we have a contradiction. □ 

Now, from this restriction, we may also find the value of x where the maximum deflection 
of u happens. To do so, we first have the following proposition. 

Proposition 2.3. Assume that A > is fixed, a G (0, 1/(1 + A)] and u(x) is a function whose 
range is contained in (0, l). Then 

2(l-u)(l-a)-\(a-u)>0 and (1 - u)(l - a) - X(a - u) > 0. (2.4a,b) 
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Proof. By collecting a terms, 

2(1 - u)(l - a) - X(a -u) = (-2(1 - u) - X)a + 2 - 2u + Au, 

which is linear in a. Therefore, since u £ (0, 1), —2(1 — u) — A < and the function takes on 
its minimum at a = 1/(1 + A). Hence, 

2(1 - U )(l - a) X(a -u)> X & ~^ Xu ) > 0. 

Similarly, one can show (1 — u)(l — a) — A(a — u) > 0. □ 

From this proposition, first integral (2.2) and Theorem 2.2, we have that the derivative of a 
solution u(- ; A, L) of (1.3) must satisfy 



u'(x; A, L) = < 



J\{\\u 


|oo-«V2(l-«)(l- 


||m||oo) - A(||u||oo - u) 




(1 - «)(1 - ||«||oo) - 


A(||m||oo - u) 


y/K\ 


liH*,- 1^2(1 -li)(l 


- ||«||oo) - A(||u||co - U) 




(l-uja-iHU) 


- A(||u||oo - u) 



x e {-L,c], 

(2.5) 

> xe(c,L), 



where ||m||oo <= (0, 1/(1 + A)] and A > 0. Next, we use derivative (2.5) to prove that the solutions 
to (1.3) must be even. 

Theorem 2.4. If u(-;X,L) is a solution of (1.3) for A > and L > 0, then u is an even 
function in x; moreover, it attains its unique maximum at x — 0, which implies that u(0) = 
||tt||oo an d u '(0) = o. 

Proof. From Lemma 2.1, we know that u attains its unique maximum at a value, say x — c, in 
(-L, L). To show that u is even, we use the fact that for every U E (0, u(c)) there exists two 
values x\ and X2 in (— L, c) and (c, L), respectively, such that u(x\) = u(x2) = U. Thus, using 
(2.5), we obtain 

(i-«)(i-NU)-A(|HU-n) du 



/Xl PL 
-L Jo 



\/X(\\u\\oo - U)y/2(1 - U)(l - ||?i||co) - A(||m||oo - u) 

(l-uja-IHUj-AdHU-u) du= r L dx 



u V A (ll u ll°o - u)^/2{\ - u)(l - IHU) ~ MpJjoo - u) 



X2 



which implies x\ + L — L — x 2 ; that is, x\ — — x 2 , proving that u is an even function in x; 
moreover, since it has only one maximum, c = 0. □ 

We summarize the results of this section in the following theorem. 
Theorem 2.5. The function u(- ;X,L) is a solution of 

- ( . U = 1 = - A < x < L; u<l, -L<x<L u'(0) = u(L) = 0, 
\^l + (u')ij (l-«) 2 ' 

(2.6) 

if and only if its even extension to (—L,L) is a solution of (1.3). Also, no solutions of (1.3) 
exist, ifu{0) > 1/(1 + A). 
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3. Time Map 

In this section, to investigate the existence of solutions of (1.3), we define and analyze the time 
map of 

I v! \ X 

x > 0; u(0) = a, u'(0) = 0, (3.1) 



\^i + (wyj (1-w) 2 ' 

foraG (0,1/(1 + A)]. 

Definition 3.1. For fixed A > 0, we define the time map T(- ; A) : (0, 1/(1 + A)] ->• R of (3.1) as 
the function that takes in a and returns the so-called "time" it takes for the solution of (3.1) to 
hit the line u — 0, i.e., it returns the smallest value L such that u(L; A, a) — 0, where u(- ; A, a) 
is the solution of (3.1). 

From this definition we have the following theorem. 

Theorem 3.2. For a given L, a certain value of X admits a solution of (2.6) if and only if we 
can find an a in (0, 1/(1 + A)] such that T(a; A) = L. 

Remark 3.3. Note that the upper bound on a comes from Theorem 2.5. 

Because of this equivalent formulation for finding solutions of (2.6), it is beneficial to have 
an analytic expression for the time map. To this end, we deduce that similar to solutions of 
(1.3), solutions of (3.1) satisfy the first integral 

1 A A 



y/1 + {u') 2 l-u 1-a 
which implies 

y/X(a - M)y/2(1 - u)(l - a) - X(a - u) 

u{x;X,a) = ry- r r , 3.2) 

(1 — u)(l — a) — A(q — u) 

for a e (0, 1/(1 + A)] and A > 0. Hence, in separating variables, integrating with respect to x 

from to L and using the change of variables u = az, we deduce from (3.2) that our time map 

is given by 

™, , \ fa f 1 (1 — az)(l — a) — Xail — z) , (' < , . , ,„ „. 

T{a;X) = J- - L . A ; v ; dz =: / K(a,z;X)dz. 3.3 

V A Jo y/l^y/2(l - az){\ -a)- Xajl - z) Jo 

Next wc prove that T(- ; A) is well defined for a e (0, 1/(1 + A)] and is in C 2 (0, 1/(1 + A)). 

Lemma 3.4. Let X > be fixed. Then T(a;X) exists for each a e (0, 1/(1 + A)]. Moreover, 
T(a; A) is differ entiable at each a G (0, 1/(1 + A)) with its derivative given by the formula 

T'(a;A)= / — (a,z;X)dz. (3.4) 
Jo da 

Furthermore, T'(a; A) is differentiable at each a G (0, 1/(1 + A)) with its derivative given by the 
formula 

T"(a;X)= f 1 ^(a,z;X)dz. (3.5) 



da 2 

Proof. First, T is well defined for a G (0, 1/(1 + A)] since 

1 

~zyj2{\ - az){l -a)- Xa{l - z) 

l + y/a 



m«; A)i < [ Vl _ : ^ 



< 



^aVTTx r 1 i a/TTa i /: 



(3.6) 
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where from inequality (2.4a) we have used the fact that 

1 



< 



^2(1 - az)(l -a)- Xa(l - z) \/Xy/T 
for a e (0, 1/(1 + A)], A > 0, z £ (0, 1). 



az 



(3.7) 



Next, 



dK 

da 



(a,z;X) 



1 



(1 — az)(l — a) — Xa(l — z) 



2y/Xa VT^ V / 2(1 - az){l -a)- Aa(l - z) 
1 + (1 - 2a)z + A(l - z) 



A y/T^Zy/2(l - az)(l -a)- Xa(l - z) 

a ((1 - az)(l - a) - Aa(l - z))(2 + A(l - z) + (2 - 4a) z) 



A 2VT^(2(l-az)(l-a)-Aa(l-z)) 3/2 
=: Hi{a,z;X) - H 2 (a, z; A) + H 3 (a, z; A), 
which is defined for a £ (0, 1/(1 + A)]. Thus, 

as: 



9a 



(a,z; A) 



< Iffxl + lfl-al + lflsl, 



so that for a e [a, 1/(1 + A)], a > 0, and z £ (0, 1), 



dK 
da 



(a,z; A) 



< 



Ci 



C 2 



V^Vl - z (1 - (1 + A)-!z) 3/2 Tl - 



(3.8) 



where Ci and C 2 are positive constants independent of z and a (see A, equations (A.l)- 
(A.4)). Therefore, we have by Theorem 10.39 in [2] that T(a; X) is differentiable at each 
a £ (0, 1/(1 + A)) with its derivative given by equation (3.4). 
Furthermore, 



2 K 



{a,z;X) 



1 



a + a\z + a 2 z 2 + a 3 z 3 



da 2 " ' ' ' ay 2 VxVT^(2(l-az)(l-a)-Xa(l-z)) 5 / 2 ' 
which is defined for each a £ (0, 1/(1 + A)) and where 

a := -1 + a(l - A), ai := a(l + 5a - 10a 2 + A + 2a 3 (2 + A)), 

a 2 := -a 3 (10 - 17a + a 2 (7 + 3A)), a 3 := a 4 (4 + 3a 2 - 2A + a(-7 + 3A)). 
Thus, for a £ [a, 1/(1 + A)], a > 0, and z £ (0, 1), 



(3.9) 



(3.10) 



8 2 K 



da 2 



{a,z;X) 



< 



1 



«i 



«2 



OS 



a 3 /VA >/T=*(2(l - az)(l - a) - Aa(l - z)) 5 /2 
(1 + A) 5 / 2 1 



< ( a + ai + a 2 + a 3 



a 3 / 2 A 3 >/T=*(l - az) 5 / 2 



< 



^/T^~z(l-{l + X)- 1 zf/ 2 ' 

where B is a positive constant that is independent of z and a. Since (1 — z) _1 / 2 (l — az)~ 5 / 2 £ 
i 1 (0, 1) for a £ (0, 1), by Theorem 10.39 in [2] we have that T'(a; A) is differentiable at each 
a £ (0, 1/(1 + A)) with its derivative given by equation (3.5). □ 

For fixed A > 0, we can now differentiate T by taking the derivative inside the integral, 
which makes it is easier to analyze the behavior of T on (0, 1/(1 + A)]. 

Proposition 3.5. Let X > be fixed. Then 
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(i) for all a e (0, 1/(1 + A)], T(a; A) > 0; 

(ii) T(a; A) as a ^ 0+; 

(iii) T'(a; A) — > +oo as a — > + ; 

(iv) lim T'(a;A)<0; 

(v) there exists a value a* <G (0, 1/(1+A)] such thatT(a*;X) = max{T(a; A) : a e (0, 1/(1 + A)]} . 
Proof. (i) This follows from definition (3.3) and Proposition 2.3. 

(ii) Since lim Q ^ + log (y^§) = 0; by inequality (3.6), we have T(a; A) — > as a — > + . 

(iii) From inequality (3.7), we have 

/ H 2 (a,z;X)dz < f \H 2 (a, z; A)| dz 
Jo Jo 

vsym r 1 i i 



< (2 + A)- 



A 2 



o Vi^Vi^ 



dz 



az 



(3.11) 



= - log^^j^O as a ^0+, 

which implies that J Q H 2 (a, z; A) dz — > as a — » + . Similarly, from inequality (A. 3), 



Jo 



i7 3 (a, z; A) dz 



< 



< 



/ \H s (a,z;X)\ 
Jo 



dz 



v^(l + A) 3 / 2 (4 + A) Z" 1 1 1 
2A 2 J 

(l + A) 3 / 2 (4 + A) ^ 



A 2 



1 - a 



o (1 - az) 3/2 y/T=. 
asa^0 + , 



: dz 



(3.12) 



and H 3 (a, z; A) dz — > as a — > + . But 



-Hi (a, z; A) dz > 



1 



2V2Aa 
1 



((1 - az)(l -a) - Xa(l - z)) dz 



l -~2 a+ 2 2 



(3.13) 



+oo as a — > + . 



2V2Aa 

Therefore, from equation (3.4) and inequalities (3.11)-(3.13), lim a _ ) . + T"(o:; A) = +oo. 
(iv) From inequality (3.8) and Theorem 10.38 of [2], we have 

A dK ( i 



lim T'(a; A) = / 



9a V 1 + A 



, z; A dz . 



(3.14) 



Now, 



dK ( 1 



da Vl + A 



,z; A = 



3A(z - 1) - (z - l) 2 - A 2 (z 2 - 2z + 3) + A 3 (z 2 + z - 1) 
XVl^XVT^(l - z + A(l + z)) 3 / 2 



Case I: If A = 1, then integrating by parts gives 



1 dK f 1 



/ da \ 1 + A 
Case II: If A ^ 1, then 

1 dK 



, z; A I dz 



o da \ 1 + A 



, z; A I dz = 



i_8 + fe-^ 8 _ 

o 4^1^ 5' 



1 + A 2 - 2A 4 + 6A 2 7A 2— n:arctan 



A(1 + A)(A-1)3 



(3.15) 



(3.16) 
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(a) If A < 1, we have VA 2 — 1 = i\J\ — A 2 and 



IX- 1 11- A i, /vTTA + \/r^A 
arctan W = arctan zW = - log 



i + A vi + A 2 b V\/tta- yr^A 

which implies from equation (3.16) that 

f 1 dK ( 1 A l + A 2 -2A 4 -3A 2 ^^^log(^g^g) 

7 9c^V1~Ta' Z j A(1 + A)(A-1)3 V "<°' fOTA<1 - 

(3.17) 

(b) If A > 1, we have < \f\ — 1 < VX+T, which implies < arctan ^^FX' anc ^ nence i 
from equation (3.16), 

- 1 dK ( 1 \ 1 + A 2 - 2A 4 + 6A 2 A /A 2_ n: arctan , 



I dK ( 1 \ l + A — zA + bA V A — 1 arctan w 

/ -TT- I z r,z;X]dz = — — — —5 < 0, for A > 1. 

J da \1 + X' J A(1 + A)(A-1) 3 

(3.18) 

Therefore, by equations (3.15), (3.17) and (3.18) the integral J* dK/da ((1 + A)" 1 , z; A) dz < 
0, for A > 0, which implies from equation (3.14) that 

lim T'(a:X)<0. 

(v) If we extend T to be defined on [0,1/(1 + A)] such that T(0; A) = 0, then by part 
(ii) T(a;X) is continuous on a compact set, [0,1/(1 + A)] and it attains it supremum in 
[0, 1/(1 + A)]. That is, by part (i), there exists a value a* e (0, 1/(1 + A)] such that T(a*; X) = 
max ae(0) i/(i + A)] T(a; A). 

□ 

A consequence of this lemma is that there exists at least one critical point of T(- ; A) in 
(0, 1/(1 + A)]. Next, we prove there exists only one. 

Proposition 3.6. Let X > be a fixed value. Then there is exactly one critical point of T in 
(0,1/(1 + A)]. 

Proof. By Proposition 3.5, we know there exists at least one critical point of T in (0,1/(1 + A)] . 
Thus, to complete the proof we only need to prove that there exists at most one critical point 
of T in (0, 1/(1 + A)]. To do so, we will show that T"(a; A) < for all a € (0, 1/(1 + A)). 

From Proposition 2.3 and equation (3.9) we see that the sign oid 2 K/da 2 for a e (0, 1/(1+A)) 
depends on the numerator 

p(z, a; A) := a (a; A) + a\(a; X)z + a 2 (a; X)z 2 + a 3 (a; X)z 3 , 

where the coefficients are defined in (3.10). Therefore, in rearranging we have 

p(z,a; X) =b Q + b\X, 

where 

6 := -1 + ha 2 z + 3a e z 3 + a(l + z) - 10a 3 (l + z)z - 7a 5 (l + z)z 2 + a A (l + Az){z + 4)z, 

61 := -a(l -£)- 3a 5 (l - z)z 2 + 2a 4 (l - z 2 )z. 

After a long computation it can be shown that bo < and b\ < for all a and z in (0, 1), 
which implies that p(z,a;X) < for all a and z in (0,1). Therefore, d 2 K/da 2 < for 
(a,z) e (0,1) x (0,1) which implies that d 2 K/da 2 < for (a,z) G (0,1/(1 + A)) x (0,1). Thus, 
by equation (3.5), T"(a; A) < for all a e (0, 1/(1 + A)). □ 
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From Proposition 3.5 and 3.6, we now have that the graph of T(a; A) looks like that of Figure 
3. Upon inspection, we see that the number of solutions of T(a; A) = L depends on two values: 
the value of T at the end point a = 1/(1 + A) and the maximum value of T for a e (0, 1/(1 + A)]. 
With this as motivation we define the functions 

M(A) := max {T(a; A) : a e (0,1/(1 + A)]} , forAe(0,oo), (3.19) 

g(X) :=T(1/(1 + A);A), forAe(0,oo), (3.20) 

and examine their properties. 




Figure 3. Plot of the time map T(a; A). 



Lemma 3.7. Let g(X) be defined as in (3.20). Then g'(X) and g"(X) exist for each A > 0. 
Moreover, g'(X) and g"(X) are given by 



f 1 86 f 1 d 2 6 

fl'(A) = y Q^(z,\)dz, and g"(X) = J ^(z,X)dz, 

z,X) :=K(l/(l + X),z;X) 



respectively, where 

A 



(i + a) 3 / 2 yr^/i - z + a(i + z) ■ 

Proof. First, 1 — z + A(l + z) = 1 + A + ( — 1 + X)z, which is linear in z, and hence, 

1 - z + A(l + z) > min{l + A,2A} > 0. (3.21) 
Then, from the definition of g(X), we have 



(3.22) 



< g(X) = — tit, \ ■ dz 

(i + A)3/2y yr^ v /i- z + A(i + z) 



f 1 i 

JO V i — 



A 

(1 + A) 3 / 2 y/mm{l + \,2\} Jo VT^z~ 
for each A > 0. Moreover, for A > 0, 

86 (A 2 -l)z + (A 2 -A + l)z 2 



8X (l + A) 5 / 2 yT^(l-z + A(l + z)) 3 / 2 

and 

_ (2A 3 - 6A - 4)z + (4A 3 - 4A 2 - 2A + 6)z 2 + (2A 3 - 4A 2 + 7A - 2)z 3 

aA 2 " (i + A) 7 /2yr^(i - z + a(i + z)) 5 /2 
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Thus, from inequality (3.21), wc obtain for each A in [a, b], where < a < 6, 

2b 2 + b + 2 1 Ci 



and 



dX 
d 2 6 



OX 2 



< 



< 



< 



(1 + a) 5 /2( m in {1 + a, 2a})3/ 2 y/T^z ~ y/\ - z 
8b 3 + 8b 2 + 156 + 12 1 _ C 2 



< 



(l + a) 7 /2( m in{l + a,2a}) 5 /2 yj\- z ~ y/T=z" 

Here C\ and Ci are positive constants independent of A and z. Since (1 — z) -1 / 2 e L 1 (0, 1), by 
Theorem 10.39 in [2] we have our result. □ 

To simplify g, we have from the definition 

A 



5(A) 



(l + A) 3/2 io a/T^vA - * + A(l + z) 



dz. 



Thus, 



9(1) 



-I 1 
4 Jo 



dz 



For A 7^ 1, the function 

A X I A ' J 

<?(A) 



(l + A) 3 / 2 7 ^((l-\)z-l) 2 -\ 2 



dz 



A 



(1 - A 2 ) 3 / 2 

A 



-8 + 1 



ds 



Vs 2 - A 2 

fl(8-l)(8*-\*)-V*ds, 



for A e (0,1), 



(1 - A 2 ) 3 / 2 | j-i J^i _ s)(A 2 - s 2 ) -1 / 2 ds, for A e (l,oo) 
/ -A(l - A 2 )- 3 / 2 (a/T^A 2 - log[(l + a/T^A 2 )/^) , for < A < 1, 
_ \A(A 2 - I)- 3 / 2 (y/W^l - sec" 1 (A)) , for 1< A < 00, 

and we have by equations (3.23)-(3.24) that 

'A(l - A 2 )" 3 / 2 (log[(l + ^1-X 2 )/X] - \/l - A 2 ) , for < A < 1, 
.9(A) = J 1/3, 

^A(A 2 - I)" 3 / 2 (a/A 



for A = 1, 
for 1 < A < 00. 



(3.23) 



(3.24) 



(3.25) 



1 — arcsec A 

Now, we may prove the following properties of g(X). 

Proposition 3.8. Let g(X) be defined as in (3.20). Then 

(i) g{X) > for all A. 

(ii) 5(A) -> as X -> 0+; 

(iii) 5(A) — > as A — > +00; 

(iv) 5' (A) +00 as A -> 0+; 

(v) g'(X) -> 0" as A -> +00; 

(vi) i/iere exisfs exactly one critical point c of g for A > 0. Moreover, c e (0, 1) and g(c) = 
L* < 00, where 

L* := max.g(A). (3.26) 

Proof. (i) Follows from inequality (3.22). 

(ii) From inequality (3.22), we see 51(A) — > as A — > + . 

(iii) From inequality (3.22), we see g(X) — > as A — > +00. 
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(iv) From equation (3.25), 



g'W = 



-(2 + A 2 )(l - A 2 + VT^V) + (1 + 2A 2 )(1 + A 2 ) log[(l + A /T^A 2 ~)/A] 



(1 - A 2 ) 5 / 2 (l + vT^A 2 ) 

for A e (0, 1), which, because log[(l+\/l — A 2 )/A] — > +oo as A — > + , implies </(A) -4- +oo 
as A -)• 0+. 
(v) From equation (3.25), 



(1 + 2A 2 )arcsec A - \/A 2 ~^T(2 + A 2 ) 1 3 
5 '(A) = i ^ 2 _ ^ '- = -—+0(X 3 ) as A ->■ +oo, 

for A > 1, which implies liniA^+co </(A) = 0. 
(vi) From equation (3.25), we deduce that g' < on (l,oo). Furthermore, by Lemma 3.7, 



,. (1 + 2A 2 )arcsec A - VW=l(2 + A 2 ) 1 
^ = * ^ (A 2 -l) 5 / 2 ^ = "is < °' 



(3.27) 



and we have g' < on [l,oo). Therefore, (/' cannot have a critical point in [l,oo). 
Next, from part (iv) of this proposition, inequality (3.27) and Lemma 3.7, we have by 
the intermediate value theorem for derivatives that there exists a c in (0, 1) such that 
g'(c) = 0. Then using equation (3.25) it can be shown that g"(X) < for A e (0,1), 
which implies that c corresponds to a local max and must be unique. Also, since g(X) > 
and g'(X) < for A g [1, oo), g(c) > g(l) > g(X) > 0, and c corresponds to a global max, 
i.e., g(c) = max A>0 .g(A). 

□ 

From this analysis we have that the graph of g(X) is given by Figure 4. Furthermore, with 
the help of Proposition 3.8 we can numerically estimate L* and find L* w 0.3499676. 



0.4 r 




Figure 4. Graph of the right end point, g(\), of the time map T(a; A). 



Next, we look at the properties of M{X). 

Proposition 3.9. Assume that M(A) is defined as in (3.19). Then 
(i) M(A) is well-defined and continuous for X > 0; 



(ii) M(A) 

(iii) M(A) 



+ as X — > +oo; 
+oo as X — >• + . 
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Proof. (i) By Proposition 3.5(v), we have that M(A) is well defined for each A > 0. Fur- 
thermore, by the definition it is easy to see that T(a; A) is continuous for each point in 
(0, 1/(1 + A)] x (0, oo ) which implies that M(A) is continuous on (0, oo). 
(ii) For any A > and a e (0, 1/(1+A)] we have by inequality (3.6) and part (i) of Proposition 
3.5 that 

0<T(a;A)< / — = dz < log' 



A 7 Vl-z ^l-il + X)-^ "A b V \/A 
which implies 

0<M(A)<^^log(l±^), for A>0. 

Now, the right hand side has the far-field behavior 2/A + 0(A~ 2 ) as A — >• +oo. Therefore, 
M(A) 0+ as A -> +oo. 
(iii) For any A > and a e (0, 1/(1 + A)] we have 

M(A) > T [ — - 1 ; A 



2(1 + A) : 



2 + 4A-z + 2A 2 (l + z) 



4\/A(l + A) 3 / 2 7o a/T^v 72 - A(-5 + i) - z + A 2 (3 + i) 

1 + 4A + 2A 2 „ , 

> — — - — = -> +oo as A H> . 

" 2VX(1 + A) 3 / 2 V2 + 5A + 4A 2 



dz 



□ 



Now, we can prove the following lemma: 
Lemma 3.10. Let T and L* be defined as (3.3) and (3.26), respectively. Then 

(i) for fixed a, T(a; A) is strictly decreasing with respect to X, which implies that M is strictly 
decreasing; 

(ii) if < L < L* , then there exist two unique constants A* and A**, where A* < A**, such 
that g(X*) = g(K„) = L; 

(iii) for each L in (0,oo), there exists a unique A* such that M(A*) = L; 

(iv) if < L < L* , then A* < A** < A*, where the values A*, A** and A* are from parts (ii) 
and (iii). 

Proof. (i) By definition (3.3), 

Tier, A) = f 1 - dz . 

V A J ^1^/2(1 - az){\ -a)- Xa{\ - z) 

For fixed values of a and z, the integrand is decreasing with respect to A, which yields 
the desired result. 

(ii) From parts (ii) and (vi) of Proposition 3.8 and < L < L* = g(a), we have by the inter- 
mediate value theorem that there exists a A* € (0, a) such that g(X*) = L. Furthermore, 
by part (vi) of Proposition 3.8, a is unique and g(a) is a maximum. This implies that g 
is increasing on (0, a), and A* is unique. 

Next, from part (iii) of Proposition 3.8, there exists an N > 0, such that X > N 
implies L > g(X), and hence, g(X) ^ L for A > N. Furthermore, since g is continuous on 
[a, N] and g(N) < L < L* , we have by the intermediate value theorem that there exists 
a A** e (&,N) such that g(A**) = L. Again using part (vi) of Proposition 3.8, we find 
that g is decreasing on (a, N] which implies that A** is unique. 

By construction, A* < a < A**, and claim (ii) holds. 
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(iii) From Proposition 3.9(H), we know there exists b > such that A > b implies M(A) < L. 
Similarly, from Proposition 3.9(iii), we know that there exists an a > (and a < b), such 
that A < a implies M(A) > L. Thus, since M is continuous on [a, b], and M(b) < L < 
M(a), by the intermediate value theorem we have that there exists an A* G [a, b] C (0, oo) 
such that M(A*) = L. 

By part (i), M(A) is strictly decreasing and thus A* is unique. 

(iv) From part (ii), A* < A**. Furthermore, M(A*) = L = g(\**) < M(A**) where the 
last inequality is from the definition of M and g — (3.19) and (3.20), respectively — and 
Proposition 3.5(iv). Thus, since M is strictly decreasing, A** < A*. 



Now, we are setup to prove the Theorem 1.1 which characterizes the solution set of (1.3). 

Proof. Let T(a; A) be defined by (3.3). Then, by Theorems 2.5 and 3.2, finding solutions to 
(1.3) is equivalent to finding a G (0, 1/(1 + A)], A > 0, such that 



Specifically, a is a solution of T(a; A) = L if and only if a — u(0;A,i), where u(- ;A,i) is a 
solution of (1.3), and the numbers of solutions of (1.3) and (3.28) are the same. Therefore, we 
look at the solutions of (3.28) in the two situations: (i) L < L*; (ii) L > L*. The following 
analysis is illustrated in Figure 5. 

(i) Let L < L* . By Lemma 3.10, we know there exists constants A*, A** and A* such that 
A* < A** < A* and g(A*) = .g(A**) = M(A*) = L. (a) If < A < A*, then A < A* which, 
since M is strictly decreasing implies M(A) > M(A*) = L. Furthermore, by the proof of 
(ii) of Lemma 3.10, g is increasing on (0, A*] and g(X) < g(A*); hence, g(X) < L < M(A), 
and by Proposition 3.5 and 3.6 there exists two solutions of T(a; A) = L. (b) If A* < 
A < A**, then by Proposition 3.8 and Lemma 3.10, g(X) G (L,g(a)]. Thus, g(X) > L and 
M(A) > L, and by Propositions 3.5 and 3.6 there exists one solution of T(a; A) — L. (c) 
If A** < A < A*, then M(A) > M(A*) = L; also by the proof of (ii) of Lemma 3.10, g 
is decreasing on [A**, A*) which implies g(X) < g(\**) = L. Hence, g(X) < L < M(A), 
and by Proposition 3.5 and 3.6 there exists two solutions of T(a: A) = L. (d) If A = A*, 
then M(\) = L, and there exists only one solution of T(a; A) = L. (c) If A > A*, then 
M (A) < M(A*) = L, and T(a; A) ^ L. 

(ii) Let L > L* which implies that g(X) < L for A > 0. Also, by Lemma 3.10, we know 
there exists a constant A* such that M(A*) = L. (a) If < A < A*, then M is strictly 
decreasing which yields M(A) > M(A*) = L. Thus, g(X) < L < M(A), and we have by 
Proposition 3.5 and 3.6 that there exists two solutions of T(a; A) = L. (b) If A = A*, 
then M(X) = L, and there exists only one solution of T(a; A) = L. (c) If A > A*, then 
M (A) < M(X*) = L, and T(a; A) ^ L. 

Now to get the first bound on A*, we first let (A,w) be a solution pair of (1.3) and consider 
the eigenvalue problem 



for ip G Hq[— L 7 L]. After multiplying (1.3) and (3.29) by <p and u, respectively, and integrating 
over (— L, L), we obtain 



□ 



T(a; A) = L. 



(3.28) 





(3.29) 
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which yields the following solvability condition for u: 



= 



X 



(1-uf 



jJLU 



dx, 



for all eigenpairs, (fj,, tp), of problem (3.29). In particular, it holds for the first eigenvalue \i\ 
of problem (3.29), which is positive and simple, and its corresponding eigenvalue ipi, which 
can — and will — be chosen to be strictly positive in [—L,L]. In choosing (fi(x) > for all x in 



[—L,L], the term A(l 



fiiu must be identically zero or change sign. It is clear that it 



is not zero; hence, there must be a value of u, say u £ (0, 1), where A(l — u)~ 2 = ^iiu. Since 
this expression must be true, we obtain < A < 4^/27. However, using Rayleigh's formula, 
we also have 



Mi 



< inf 




dx:<p€Hi[-L,L], |M| 2 = 1 



dx:ipeH^[-L,L], \\tph = l 



where K\ = 7r 2 /(4L 2 ) is the first eigenvalue of —d 2 /dx 2 on [-L, L] with homogeneous Dirichlet 
boundary conditions and || • || 2 is the standard L 2 norm on \—L,L\. Therefore, the value 
A < 4^i/27 < 4ki/27 = n 2 /(27L 2 ). 

To get the second bound we note that if u is a solution to (1.3), then (1 — u)~ 2 > 1. Therefore, 
integrating the differential equation in (1.3) from to a; and using the aforementioned inequality 
gives 

u'(x) 



^l + u'{xf 



> Xx, 



< x < L, 



where we have used that fact that u'(0) = 0. From Lemma 2.1, we know that — u'(x) — \u'(x)\; 
thus, 



Xx < 



\u'(x)\ 



which upon taking x — > L~ yields A < L^ 1 . 
Hence, A* < min{L-\ n 2 /(27L 2 )}. 



< 1, 



□ 




(a) L < L* 




Figure 5. Graphical illustration of Theorem 1.1 
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4. Conclusion 

In this work, we have analyzed the solution set of the one-dimensional prescribed mean curva- 
ture problem 

/ \ ' 

1 "' '' -L<x< L; u{-L) = u(L) = 0, (1.3) 



with u < 1 in [—L,L], for positive L and positive A. In particular, we have shown that 
the solution set undergoes two bifurcations: a saddle node bifurcation and codimension 2 
splitting bifurcation, which depends on the parameters L and A. As a result of this analysis, 
we speculate that the latter, which is not present in the corresponding semilinear problem 
(—Ait = A _1 (l — u)~ 2 with the same boundary conditions), is due to the interplay of the mean 
curvature operator and finite-singularity forcing. This is because while other one-dimensional 
prescribed mean curvature equations have exhibited a disappearing solutions behavior (c.f. 
[57]), none with a continuous forcing have exhibited this splitting bifurcation. With this in mind 
it would be desirable to fully characterize the solution set of the one-dimensional prescribed 
mean curvature equation with f(u) = A(l — u)~ p , for p > 0, and see if the splitting bifurcation 
generalizes to this class of problems. 

Appendix. Bounding Hi, for i = 1,2,3. 

By Proposition 2.3 for a e [a, 1/(1 + A)], a > 0, and z e (0, 1), we deduce 

1 (1 — az)(\ — a) — Aa(l — z) 



< Hi(a, z; A) = 



< 



2v / Aa yi~^^2(l - az)(l -a)- \a(l - z) 
1 1 



0< \H 2 {a,z;X)\ 



2y / Aa 71^^/2(1 - az)(l -a)- \a(l - z) ' 
and 

|l + (l-2a)z + A(l-z)| 
A 71^^2(1 - az)(l -a)- Xa(l - z) 
fa l + |l-2a|+A 

A 71^^2(1 - az)(l -a)- Xaji - z) 
2 + A 

A 71^^2(1 - az)(l - a) - Xa{l - z) ' 
Moreover from Proposition 2.3, 

A A^cy 
2(1 - az)(l -a)- Xa(l - z) > T^C 1 - az + Xaz) > -^~\ z > °' 



< 



< 



which implies 



and 



v/m i vm i 

- 1 < 2A3/2 a ^VT^~z < 2A3/2 a y/zVT^t ( ' j 



0< |g2 l<(2 + A) ^^^_ . (A.2) 
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Again, from Proposition 2.3 for a e [a, 1/(1 + A)], a > 0, and z e (0, 1), we have 

fa((l - az)(l -a)- Aa(l - z))\2 + A(l -z) + {2- Aa)z\ 



Q<\Hz(a,z;\)\ 



2>/r=z(2(l - az){\ -a)- Aa(l - z)f /2 
2 + A + 2|l-2a| 



< 



Also, from Proposition 2.3, 



* 2^1^(2(1 - az)(l -a)- Aa(l - z)) 3/2 
4 + A 

A 2^1^ (2(1 - az)(l -a)- Aa(l - z)) 3/2 ' 



2(1 - az)(l - a) - Aa(l - z) > -(1 - az + Xaz) > (1 - az) > 0, 

1 -|- A 1 — |— A 



-^-r(l — az + Xaz) > 
+ A V ' ~ 1 + 

which implies 

V^(l + A) 3 /2(4 + A) 1 1 



0< \Ho\ < 



< 



2A 2 (i _ az f 2 VT^ 

(1 + A)(4 + A) 1 1 



2A 2 (l-(l + X)^zf /2 VT^' 
Furthermore, for A > 

f 1 I f 1 1 1 

/ F n dz = n > / a. dz^2(l + A- 1 ), 

7o a/^a/T^ 7o (1 - (1 + A)-!z) 3/2 VT^ V ; 

which means that l/(vWl - z) and (1 - (1 + X)~ 1 z)~ 3 / 2 /y/l - z arc in L x (0, 1). 



(A.3) 



(A.4) 
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